Phenotyping the virulence of SARS-CoV-2 variants in hamsters by digital pathology and machine learning (preprint)

SARS-CoV-2 has continued to evolve throughout the COVID-19 pandemic, giving rise to multiple variants of concern (VOCs) with different biological properties. As the pandemic progresses, it will be essential to test in near real time the potential of any new emerging variant to cause severe disease. BA.1 (Omicron) was shown to be attenuated compared to the previous VOCs like Delta, but it is possible that newly emerging variants may regain a virulent phenotype. Hamsters have been proven to be an exceedingly good model for SARS-CoV-2 pathogenesis. Here, we aimed to develop robust quantitative pipelines to assess the virulence of SARS-CoV-2 variants in hamsters. We used various approaches including RNAseq, RNA in situ hybridization, immunohistochemistry, and digital pathology, including software assisted whole section imaging and downstream automatic analyses enhanced by machine learning, to develop methods to assess and quantify virus-induced pulmonary lesions in an unbiased manner. Initially, we used Delta and Omicron to develop our experimental pipelines. We then assessed the virulence of recent Omicron sub-lineages including BA.5, XBB, BQ.1.18, BA.2 and BA.2.75. We show that in experimentally infected hamsters, accurate quantification of alveolar epithelial hyperplasia and macrophage infiltrates represent robust markers for assessing the extent of virus-induced pulmonary pathology, and hence virus virulence. In addition, using these pipelines, we could reveal how some Omicron sub-lineages (e.g., BA.2.75) have regained virulence compared to the original BA.1. Finally, to maximise the utility of the digital pathology pipelines reported in our study, we developed an online repository containing representative whole organ histopathology sections that can be visualised at variable magnifications (https://covid-atlas.cvr.gla.ac.uk). Overall, this pipeline can provide unbiased and invaluable data for rapidly assessing newly emerging variants and their potential to cause severe disease.

Inconsistent directions of change in case severity across successive SARS-CoV-2 variant waves suggests an unpredictable future (preprint)

ObjectiveTo determine how the severity of successively dominant SARS-CoV-2 variants changed over the course of the COVID-19 pandemic. DesignRetrospective cohort analysis. SettingCommunity- and hospital-sequenced COVID-19 cases in the NHS Greater Glasgow and Clyde (NHS GG&C) Health Board. ParticipantsAll sequenced non-nosocomial adult COVID-19 cases in NHS GG&C infected with the relevant SARS-CoV-2 lineages during analysis periods. B.1.177/Alpha: 1st November 2020 - 30th January 2021 (n = 1640). Alpha/Delta: 1st April - 30th June 2021 (n = 5552). AY.4.2 Delta/non-AY.4.2 Delta: 1st July - 31st October 2021 (n = 9613). Non-AY.4.2 Delta/Omicron: 1st - 31st December 2021 (n = 3858). Main outcome measuresAdmission to hospital, ICU, or death within 28 days of positive COVID-19 test ResultsFor B.1.177/Alpha, 300 of 807 B.1.177 cases were recorded as hospitalised or worse, compared to 232 of 833 Alpha cases. After adjustment, the cumulative odds ratio was 1.51 (95% CI: 1.08-2.11) for Alpha versus B.1.177. For Alpha/Delta, 113 of 2104 Alpha cases were recorded as hospitalised or worse, compared to 230 of 3448 Delta cases. After adjustment, the cumulative odds ratio was 2.09 (95% CI: 1.42-3.08) for Delta versus Alpha. For non-AY.4.2 Delta/AY.4.2 Delta, 845 of 8644 non-AY.4.2 Delta cases were recorded as hospitalised or worse, compared to 101 of 969 AY.4.2 Delta cases. After adjustment, the cumulative odds ratio was 0.99 (95% CI: 0.76-1.27) for AY.4.2 Delta versus non-AY.4.2 Delta. For non-AY.4.2 Delta/Omicron, 30 of 1164 non-AY.4.2 Delta cases were recorded as hospitalised or worse, compared to 26 of 2694 Omicron cases. After adjustment, the median cumulative odds ratio was 0.49 (95% CI: 0.22-1.06) for Omicron versus non-AY.4.2 Delta. ConclusionsThe direction of change in disease severity between successively emerging SARS-CoV-2 variants of concern was inconsistent. This heterogeneity demonstrates that severity associated with future SARS-CoV-2 variants is unpredictable.

The hyper-transmissible SARS-CoV-2 Omicron variant exhibits significant antigenic change, vaccine escape and a switch in cell entry mechanism (preprint)

Vaccination-based exposure to spike protein derived from early SARS-CoV-2 sequences is the key public health strategy against COVID-19. Successive waves of SARS-CoV-2 infections have been characterised by the evolution of highly mutated variants that are more transmissible and that partially evade the adaptive immune response. Omicron is the fifth of these Variants of Concern (VOCs) and is characterised by a step change in transmission capability, suggesting significant antigenic and biological change. It is characterised by 45 amino acid substitutions, including 30 changes in the spike protein relative to one of the earliest sequences, Wuhan-Hu-1, of which 15 occur in the receptor-binding domain, an area strongly associated with humoral immune evasion. In this study, we demonstrate both markedly decreased neutralisation in serology assays and real-world vaccine effectiveness in recipients of two doses of vaccine, with efficacy partially recovered by a third mRNA booster dose. We also show that immunity from natural infection (without vaccination) is more protective than two doses of vaccine but inferior to three doses. Finally, we demonstrate fundamental changes in the Omicron entry process in vitro, towards TMPRSS2-independent fusion, representing a major shift in the replication properties of SARS-CoV-2. Overall, these findings underlie rapid global transmission and may alter the clinical severity of disease associated with the Omicron variant.

The SARS-CoV-2 Alpha variant is associated with increased clinical severity of disease (preprint)

BackgroundThe B.1.1.7 (Alpha) SARS-CoV-2 variant of concern was associated with increased transmission relative to other variants present at the time of its emergence and several studies have shown an association between the B.1.1.7 lineage infection and increased 28-day mortality. However, to date none have addressed the impact of infection on severity of illness or the need for oxygen or ventilation. MethodsIn this prospective clinical cohort sub-study of the COG-UK consortium, 1475 samples from hospitalised and community cases collected between the 1st November 2020 and 30th January 2021 were collected. These samples were sequenced in local laboratories and analysed for the presence of B.1.1.7-defining mutations. We prospectively matched sequence data to clinical outcomes as the lineage became dominant in Scotland and modelled the association between B.1.1.7 infection and severe disease using a 4-point scale of maximum severity by 28 days: 1. no support, 2. oxygen, 3. ventilation and 4. death. Additionally, we calculated an estimate of the growth rate of B.1.1.7-associated infections following introduction into Scotland using phylogenetic data. ResultsB.1.1.7 was responsible for a third wave of SARS-CoV-2 in Scotland, and rapidly replaced the previously dominant second wave lineage B.1.177) due to a significantly higher transmission rate ([~]5 fold). Of 1475 patients, 364 were infected with B.1.1.7, 1030 with B.1.177 and 81 with other lineages. Our cumulative generalised linear mixed model analyses found evidence (cumulative odds ratio: 1.40, 95% CI: 1.02, 1.93) of a positive association between increased clinical severity and lineage (B.1.1.7 versus non-B.1.1.7). Viral load was higher in B.1.1.7 samples than in non-B.1.1.7 samples as measured by cycle threshold (Ct) value (mean Ct change: -2.46, 95% CI: -4.22, -0.70). ConclusionsThe B.1.1.7 lineage was associated with more severe clinical disease in Scottish patients than co-circulating lineages. FundingCOG-UK is supported by funding from the Medical Research Council (MRC) part of UK Research & Innovation (UKRI), the National Institute of Health Research (NIHR) and Genome Research Limited, operating as the Wellcome Sanger Institute. Funding was also provided by UKRI through the JUNIPER consortium (grant number MR/V038613/1). Sequencing and bioinformatics support was funded by the Medical Research Council (MRC) core award (MC UU 1201412).

SARS-CoV-2 lineage B.1.1.7 is associated with greater disease severity among hospitalised women but not men (preprint)

Background Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineage B.1.1.7 has been associated with an increased rate of transmission and disease severity among subjects testing positive in the community. Its impact on hospitalised patients is less well documented. Methods We collected viral sequences and clinical data of patients admitted with SARS-CoV-2 and hospital-onset COVID-19 infections (HOCIs), sampled 16/11/2020 - 10/01/2021, from eight hospitals participating in the COG-UK-HOCI study. Associations between the variant and the outcomes of all-cause mortality and intensive therapy unit (ITU) admission were evaluated using mixed effects Cox models adjusted by age, sex, comorbidities, care home residence, pregnancy and ethnicity. Results Sequences were obtained from 2341 inpatients (HOCI cases = 786) and analysis of clinical outcomes was carried out in 2147 inpatients with all data available. The hazard ratio (HR) for mortality of B.1.1.7 compared to other lineages was 1.01 (95% CI 0.79-1.28, P=0.94) and for ITU admission was 1.01 (95% CI 0.75-1.37, P=0.96). Analysis of sex-specific effects of B.1.1.7 identified increased risk of mortality (HR 1.30, 95% CI 0.95-1.78) and ITU admission (HR 1.82, 95% CI 1.15-2.90) in females infected with the variant but not males (mortality HR 0.82, 95% CI 0.61-1.10; ITU HR 0.74, 95% CI 0.52-1.04). Conclusions In common with smaller studies of patients hospitalised with SARS-CoV-2 we did not find an overall increase in mortality or ITU admission associated with B.1.1.7 compared to other lineages. However, women with B.1.1.7 may be at an increased risk of admission to intensive care and at modestly increased risk of mortality.

Genetic epidemiology of SARS-CoV-2 transmission in renal dialysis units - a high risk community-hospital interface (preprint)

Objectives: Patients requiring haemodialysis are at increased risk of serious illness with SARS-CoV-2 infection. To improve the understanding of transmission risks in six Scottish renal dialysis units, we utilised the rapid whole-genome sequencing data generated by the COG-UK consortium. Methods: We combined geographical, temporal and genomic sequence data from the community and hospital to estimate the probability of infection originating from within the dialysis unit, the hospital or the community using Bayesian statistical modelling and compared these results to the details of epidemiological investigations. Results: Of 671 patients, 60 (8.9%) became infected with SARS-CoV-2, of whom 16 (27%) died. Within-unit and community transmission were both evident and an instance of transmission from the wider hospital setting was also demonstrated. Conclusions: Near-real-time SARS-CoV-2 sequencing data can facilitate tailored infection prevention and control measures, which can be targeted at reducing risk in these settings. Key words: SARS-CoV-2, COVID-19, haemodialysis, renal dialysis unit, infection control, rapid sequencing, outbreak, nosocomial Key words: SARS-CoV-2, COVID-19, haemodialysis, renal dialysis unit, infection control, rapid sequencing, outbreak, nosocomial

In vitro evolution of Remdesivir resistance reveals genome plasticity of SARS-CoV-2 (preprint)

Remdesivir (RDV) is used widely for COVID-19 patients despite varying results in recent clinical trials. Here, we show how serially passaging SARS-CoV-2 in vitro in the presence of RDV selected for drug-resistant viral populations. We determined that the E802D mutation in the RNA-dependent RNA polymerase was sufficient to confer decreased RDV sensitivity without affecting viral fitness. Analysis of more than 200,000 sequences of globally circulating SARS-CoV-2 variants show no evidence of widespread transmission of RDV-resistant mutants. Surprisingly, we also observed changes in the Spike (i.e., H69 E484, N501, H655) corresponding to mutations identified in emerging SARS-CoV-2 variants indicating that they can arise in vitro in the absence of immune selection. This study illustrates SARS-CoV-2 genome plasticity and offers new perspectives on surveillance of viral variants.

Epidemic waves of COVID-19 in Scotland: a genomic perspective on the impact of the introduction and relaxation of lockdown on SARS-CoV-2 (preprint)

The second SARS virus, SARS-CoV-2, emerged in December 2019, and within a month was globally distributed. It was first introduced into Scotland in February 2020 associated with returning travellers and visitors. By March it was circulating in communities across the UK, and to control COVID-19 cases, and prevent overwhelming of the National Health Service (NHS), a 'lockdown' was introduced on 23rd March 2020 with a restriction of people's movements. To augment the public health efforts a large-scale genome epidemiology effort (as part of the COVID-19 Genomics UK (COG-UK) consortium) resulted in the sequencing of over 5000 SARS-CoV-2 genomes by 18th August 2020 from Scottish cases, about a quarter of the estimated number of cases at that time. Here we quantify the geographical origins of the first wave introductions into Scotland from abroad and other UK regions, the spread of these SARS-CoV-2 lineages to different regions within Scotland (defined at the level of NHS Health Board) and the effect of lockdown on virus 'success'. We estimate that approximately 300 introductions seeded lineages in Scotland, with around 25% of these lineages composed of more than five viruses, but by June circulating lineages were reduced to low levels, in line with low numbers of recorded positive cases. Lockdown was, thus, associated with a dramatic reduction in infection numbers and the extinguishing of most virus lineages. Unfortunately since the summer cases have been rising in Scotland in a second wave, with >1000 people testing positive on a daily basis, and hospitalisation of COVID-19 cases on the rise again. Examining the available Scottish genome data from the second wave, and comparing it to the first wave, we find that while some UK lineages have persisted through the summer, the majority of lineages responsible for the second wave are new introductions from outside of Scotland and many from outside of the UK. This indicates that, while lockdown in Scotland is directly linked with the first wave case numbers being brought under control, travel-associated imports (mostly from Europe or other parts of the UK) following the easing of lockdown are responsible for seeding the current epidemic population. This demonstrates that the impact of stringent public health measures can be compromised if following this, movements from regions of high to low prevalence are not minimised.

The circulating SARS-CoV-2 spike variant N439K maintains fitness while evading antibody-mediated immunity (preprint)

SARS-CoV-2 can mutate to evade immunity, with consequences for the efficacy of emerging vaccines and antibody therapeutics. Herein we demonstrate that the immunodominant SARS-CoV-2 spike (S) receptor binding motif (RBM) is the most divergent region of S, and provide epidemiological, clinical, and molecular characterization of a prevalent RBM variant, N439K. We demonstrate that N439K S protein has enhanced binding affinity to the hACE2 receptor, and that N439K virus has similar clinical outcomes and in vitro replication fitness as compared to wild- type. We observed that the N439K mutation resulted in immune escape from a panel of neutralizing monoclonal antibodies, including one in clinical trials, as well as from polyclonal sera from a sizeable fraction of persons recovered from infection. Immune evasion mutations that maintain virulence and fitness such as N439K can emerge within SARS-CoV-2 S, highlighting the need for ongoing molecular surveillance to guide development and usage of vaccines and therapeutics.

Evaluating the effects of SARS-CoV-2 Spike mutation D614G on transmissibility and pathogenicity (preprint)

Global dispersal and increasing frequency of the SARS-CoV-2 Spike protein variant D614G are suggestive of a selective advantage but may also be due to a random founder effect. We investigate the hypothesis for positive selection of Spike D614G in the United Kingdom using more than 25,000 whole genome SARS-CoV-2 sequences. Despite the availability of a large data set, well represented by both Spike 614 variants, not all approaches showed a conclusive signal of positive selection. Population genetic analysis indicates that 614G increases in frequency relative to 614D in a manner consistent with a selective advantage. We do not find any indication that patients infected with the Spike 614G variant have higher COVID-19 mortality or clinical severity, but 614G is associated with higher viral load and younger age of patients. Significant differences in growth and size of 614G phylogenetic clusters indicate a need for continued study of this variant.

Genomic epidemiology of SARS-CoV-2 spread in Scotland highlights the role of European travel in COVID-19 emergence (preprint)

SARS-CoV-2, the causative agent of COVID-19, emerged in Wuhan, China in December 2019 and spread rapidly throughout the world. Understanding the introductions of this new coronavirus in different settings may assist control efforts and the establishment of frameworks to support rapid response in future infectious disease outbreaks. We investigated the first four weeks of emergence of the SARS-CoV-2 virus in Scotland after the first case reported on the 1st March 2020. We obtained full genome sequences from 452 individuals with a laboratory-confirmed diagnosis of COVID-19, representing 20% of all cases until 1st April 2020 (n=2310). This permitted a genomic epidemiology approach to study the introductions and spread of the SARS-2 virus in Scotland. From combined phylogenetic and epidemiological analysis, we estimated at least 113 introductions of SARS-CoV-2 into Scotland during this period. Clusters containing multiple sequences suggestive of onward transmission occurred in 48/86 (56%). 42/86 (51%) clusters had no known international travel history indicating undetected introductions. The majority of viral sequences were most closely related to those circulating in other European countries, including Italy, Austria and Spain. Travel-associated introductions of SARS-CoV-2 into Scotland predated travel restrictions in the UK and other European countries. The first local transmission occurred three days after the first case. A shift from travel-associated to sustained community transmission was apparent after only 11 days. Undetected introductions occurred prior to the first known case of COVID-19. Earlier travel restrictions and quarantine measures might have resulted in fewer introductions into Scotland, thereby reducing the number of cases and the subsequent burden on health services. The high number of introductions and transmission rates were likely to have impacted on national contact tracing efforts. Our results also demonstrate that local real-time genomic epidemiology can be used to monitor transmission clusters and facilitate control efforts to restrict the spread of COVID-19. FundingMRC (MC UU 1201412), UKRI/Wellcome (COG-UK), Wellcome Trust Collaborator Award (206298/Z/17/Z - ARTIC Network; TCW Wellcome Trust Award 204802/Z/16/Z Research in contextO_ST_ABSEvidence before this studyC_ST_ABSCoronavirus disease-2019 (COVID-19) was first diagnosed in Scotland on the 1st of March 2020 following the emergence of the causative severe acute respiratory system coronavirus 2 (SARS-CoV-2) virus in China in December 2019. During the first month of the outbreak in Scotland, 2310 positive cases of COVID-19 were detected, associated with 1832 hospital admissions, 207 intensive care admissions and 126 deaths. The number of introductions into Scotland and the source of those introductions was not known prior to this study. Added value of this studyUsing a combined phylogenetic and epidemiological approach following real-time next generation sequencing of 452 SARS-CoV-2 samples, it was estimated that the virus was introduced to Scotland on at least 113 occasions, mostly from other European countries, including Italy, Austria and Spain. Localised outbreaks occurred in the community across multiple Scottish health boards, within healthcare facilities and an international conference and community transmission was established rapidly, before local and international lockdown measures were introduced.